The \JTfit Collaboration Average of D meson mixing data: Spring 2012 

(UTfit Collaboration^ 



O 



Oh 



> 

O 
(N 



X 



A.J. Bevan/ M. Bona/ M. Ciuchini,^ D. Derkach,^ E. Franco,^ V. Lubicz,^'^ G. Martinelli,'*^ ^ F. Parodi/ 
M. Pierini,^ C. SchiaviJ L. Silvestrini,'^ V. Sordini,® A. Stocchi,^ C. Tarantino,^' ^ and V. Vagnoni^° 

^ Queen Mary, University of London, Mile End Road, London El iNS, United Kingdom 
^INFN, Sezione di Roma Tre, Via della Vasca Navale 84, LOOI46 Roma, Haly 
^CERN, CH-1211 Geneva 23, Switzerland 
'^INFN, Sezione di Roma, Piazzale A. Mora 2, 1-00185 Roma, Italy 
^ Dipartimento di Fisica, Universita di Roma Tre, 
Via della Vasca Navale 84, FOOI46 Roma, Italy 
^SISSA-ISAS, Via Bonomea 265, 1-34136 Trieste, Italy 
Dipartimento di Fisica, Universita di Genova and INFN, Via Dodecaneso 33, I-I6I46 Genova, Italy 
^IPNL-IN2P3, 4 Rue Enrico Fermi, F-69622 Villeurhanne Cedex, France 
^ Laboratoire de VAccelerateur Lineaire, IN2P3-CNRS et Universite de Pans-Sud, BP 34, F-91898 Orsay Gedex, France 

INFN, Sezione di Bologna, Via Irnerio 46, 1-40126 Bologna, Italy 

We derive constraints on the parameters M12, ri2 and $12 that describe D meson mixing using 
all available data, allowing for CP violation. We also provide posterior distributions and predictions 
for observable parameters appearing in D physics. 

Meson-antimeson mixing in the neutral D system has been established only in 2007 [TH3]- Early combinations of 
available data allowed to put stringent constraints on New Physics (NP) contributions, although the possibility of 
non-standard CP violation remained open More recently, CP violation in the D system received considerable 

attention after the measurement at hadron colliders of large direct CP violation in — >■ tttt and D — )■ KK decays 
[HI [TU] , which may signal the presence of NP [TT| - [TB] . It then becomes crucial to extract updated information on 
the mixing amplitude in order both to disentangle more precisely indirect and direct CP violation in D — >■ tttt and 
D — >■ KK, and to obtain up-to-date constraints on NP in AC — 2 transitions that can be used to constrain NP 
contributions to AC — 1 processes in any given model. 

In this letter, we perform a fit to the experimental data in Table |T] following the statistical method described in 
ref. ISnj- We assume that all Cabibbo allowed (and doubly Cabibbo suppressed) decay amplitudes in the phase 
convention CPlZ?) = \D) and CP|/) — ?7cp|/) satisfy the relation A{D f) = ri(^pA{D /), which is expected to 
hold in the SM (in the standard CKM phase convention) with an accuracy much better than present experimental 
errors. In the same approximation this implies real. For singly Cabibbo suppressed decays D*^ K^ K^ and 



we allow for direct CP violation to be present. 



and \q/p\, with \Dl,s) = p\D") ± q\D°) and |pp 
terms of x, y and \q/p\ using the following formulae [H 

' = arg(g/p) = arg 



We assume flat priors for x = Amo/^D, U = ATd/C^^^d) 
1. We can then express all mixing-related observables in 
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Observable 


Value 
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0.00978 


m 


2\fRr) cos (Sktt 
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0.0035 1 


0.0471 


m 


1\JRd sin (J/t'ir 
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1 


m 




(0.3030 ±0.0189)% 


1 


0.77 


-0.87 




m 




(-0.024 ± 0.052)% 


0.77 


1 


-0.94 




m 


(y+)if7r 


(0.98 ±0.78)% 


-0.87 


-0.94 


1 




m 


Ad 


(-2.1 ± 5.4)% 


1 


0.77 


-0.87 




m 
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CP asymmetry 


Value 






A(i)/rco 




Reference 


Acp{D° K+K- 


) (-0.24 ±0.24)% 










EiEzl 


Acp{D" ^ TT+TV-) 


(0.11 ±0.39)% 
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(-0.82 ±0.21 ±0.11)% 
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TABLE I. Experimental data used in the analysis of D mixing, from ref. |38) . a = (1 ± \q/p\)^/2 and AAcp ~ Acp{D'^ — >■ 
K~) — Acp{D'^ — > tt^tt"). Asymmetric errors have been symmetrized. We do not use measurements that do not allow for 
CP violation in mixing, except for ref. 27 (as shown in ref. the results for x and y from the Dalitz analysis of D — >■ Kann 
are not sensitive to the assumptions about CP violation in mixing). 



with Sf a strong phase and Aj^ forced to vanish in the fit. In addition, for the CP asymmetries we have 

where Df{t) is the observed distribution of proper decay time and tjjo is the lifetime of the neutral D mesons. 

For the purpose of constraining NP, it is useful to express the fit results in terms of the AC = 2 effective Hamiltonian 
matrix elements M12 and ri2: 



1 jx^ + S^if + ^ _ |ri2|2+4|Afi2|2-(a;2±y2)|q/p|2/^2 



T£,\/4(l-(5^) TD V 1-6^ 4|Mi2ri2| 
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parameter 


result @ 68% prob. 


95% prob. range 


IA/12I [1/ps] 


(6.9 ±2.4) ■ 10"^ 


[2.1,11.5] ■ 10"^ 


|ri2i [1/ps] 


(17.2 ±2.5) • 10"^ 


[12.3,22.4] • 10"^ 


$12 n 


(-6 ±9) 


[-37, 13] 


X 


(5.6 ±2.0) • 10"=* 


[1.4,9.6] • 10"^ 
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(7.0 ± 1.0) • 10-=* 


[5.0,9.1] ■ 10"^ 




(5.3 ±7.7) • 10^2 


[-8.5,25.6] • 10^2 


<^n 


(-2.4 ± 2.9) 


[-8.8,3.7] 
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(0.7 ±0.8) • 10-=* 


[-0.9,2.3] • 10"=^ 
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(11 ± 14) • 10"^ 


[-15,44] • 10"^ 


Rm 


(4.0 ± 1.4) ■ 10"^ 


[1.7, 7.2] ■ 10"^ 
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(3.27 ±0.08) • 10"=^ 


[3.10,3.44] • 10"^ 


Sktt [°] 


(18± 12) 


[-14,40] 


SktvtvO [°] 


(31 ± 20) 


[-11,73] 




) (-2.6 ±2.2) • 10~^ 


[-7.1,1.9] ■ 10"^ 




(4.1 ±2.4) • 10"^ 


[-0.8,9.0] ■ 10"=^ 




(6.6 ± 1.6) • 10"^ 


[-9.8,3.5] • 10"^ 



TABLE II. Results of the fit to D mixing data. Aa^^ = acp(-D" ^ K+K-) - a'^^(D° tt+tt"). 



with $12 = argri2/Mi2. Consistently with the assumption A{D —>■/)— A{D — > /), ri2 can be taken real with 
negligible NP contributions, and a nonvanishing $12 can be interpreted as a signal of new sources of CP violation in 
M12. For the sake of completeness, we report here also the formulae to compute the observables x, y and S from M12 
and ri2: 



V2Am = sign(cos$i2)\/4|Mi2|2 - |ri2|2 ± V(4|Afi2|2 + |ri2|2)2 - 16|Afi2|2|ri2|2sin2$ 



V2 AF - y |Fi2|2 - 4|Mi2|2 + Y/(4|Afi2|2 ± |Fi2|2)2 - 16|Mi2|2|Fi2|2 sin^ Cf>i2 , 

2|Mi2||Fi2|sin$i2 

(Ato)2 + |Fi2|2 ' 

in agreement with [42| up to a factor of V2. 
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FIG. 1. One-dimensional p.d.f. for the parameters |Afi2|, |ri2| and <l>i 



The results of the fit are reported in Tabl e [III The corresponding p.d.f are shown in Figs. [T] and [2] Some two- 
dimensional correlations are displayed in Fig. I3j 

A direct comparison with the HFAG results [38] is not straightforward, as our fit docs not fall into any of the 
HFAG categories (no CPV, no direct CPV, direct CPV) , since we allow for direct CP violation only in singly Cabibbo 
suppressed decays. However, our fit results should be close to the "no direct CPV" HFAG fit. Indeed, we find 
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MJps'] MJps'] X |q/p|-1 



FIG. 3. Two-dimensional p.d.f. for |ri2| vs |A'fi2| (top left), "l>i2 vs |Afi2l (top right), y ys x (bottom left) and vs \q/p\ — 1 
(bottom right). 



compatible results within errors. We notice, however, that HFAG performs a fit with four independent parameters 
{x, y, (j> and |(7/p|), while only three of these parameters are independent, as can be seen from eq. M. In particular, 
(f> should vanish for \q/p\ = 1. This feature can be seen in Fig. [s] (up to the smoothing of the p.d.f) but not in the 
equivalent plot from HFAG, which displays completely different 2-dimensional contours. We can but recommend that 
in the future HFAG takes the relation (p — aig{y + iSx) always into account. 

The resuhs in Table [h] can be used to constrain NP contributions to D ~ D mixing and decays. 
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